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We propose a method for calculating DIS jet production cross sections in QCD at NLO accuracy 
with consistent treatment of heavy quarks. The scheme relies on the dipole subtraction method 
for jets, which we extend to all possible initial state splittings with heavy partons, so that the 
Aivazis-Collins-Olness-Tung massive coUinear factorization scheme (ACOT) can be applied. As a 
first check of the formalism we recover the ACOT result for the heavy quark structure function 
using a dedicated Monte Carlo program. 



I. INTRODUCTION 

The contributions of heavy quarks to the inclusive DIS 
processes can be consistently described using so called 
general mass schemes (GM) (see e.g. pL^ for a review). 
The GM schemes are designed to solve the problem oc- 
curring in perturbative QCD when two or more large 
scales change significantly with respect to each other. In 
DIS the large scales are the virtuality of exchanged bo- 
son, and the mass toq of the heavy quark Q (assume 
for simplicity that there is only one heavy quark). Let 
us thus consider the ratio p = rriQ/Q^. There are two 
distinct regions, p ~ and p ^ 1, where different ap- 
proaches are used. In the first region Q is treated as 
massless and the factorization theorem is used to resum 
the resulting coUinear singularities into the parton distri- 
bution function (PDF) for Q. In the other region {p ^ 1) 
the effects of the heavy quark mass are fully taken into 
account, order by order in the perturbation expansion. 
In particular, the arising powers of as log p are finite and 
no resummation is required. This approach is referred to 
as 'fixed order'. 

The GM schemes provide a description over the whole 
range of p, by means of different resummation methods of 
the potentially large {as log p)" terms. There are several 
solutions of this kind [2 - 5] which are used in the descrip- 
tion of inclusive data. For the purpose of our investi- 
gation the ACOT scheme |2] is of particular interest as 
it actually is a factorization theorem with heavy quarks, 
proved to all orders in QCD [6 . 

The GM schemes have been so far formulated and used 
for inclusive DIS processes only. The purpose of this 
letter is to report on a new ACOT-based GM scheme 
for jet production in DIS. All necessary calculations were 
made in [T and a detailed analysis will be presented in a 
forthcoming paper. 

Analysis of the jet production processes is a key tool to 
study various aspects of perturbative QCD. On the theo- 
retical side the calculations are much more involved than 
in the inclusive case. This is due to infrared (IR) singu- 
larities that appear at intermediate stages of calculations 



but eventually get canceled in the physical cross sections. 
The problems originate in the fact that the phase space 
is not complete as the final state partons have to pro- 
duce distinct jets. Thus, in practice, only Monte Carlo 
(MC) methods are applicable here and special methods 
for canceling IR singularities between real and virtual 
corrections have to be used. 

One of the exact solutions is the dipole subtraction 
method (DSM) [HI E]- It was constructed initially for 
massless quarks only, but later a complete extension for 
massive quarks in the final state was made |10l . The 
method made use of Ref. [12j, where photon radiation 
off (possibly massive) fermions was considered, also in 
the initial state. 

It is important to note that the GM schemes deal with 
initial state splitting processes, where two of the partic- 
ipating partons are heavy quarks. None of the methods 
mentioned above take this into account completely. In 
Ref. [TT] there are no massive partons in the initial state 
splittings, while Ref. [T^ considers only Q — ^ Q7 emis- 
sions. 

The method presented in this paper consists in two ba- 
sic elements: i) dipole subtraction method extended to 
all possible QCD splitting processes with heavy quarks, 
ii) the ACOT factorization scheme for initial state quasi- 
collinear singularities. The last notion corresponds to a 
situation, where both transverse momentum and mass of 
the emitted parton tend to zero in a uniform way [13' . In 
what follows we shall refer to these quasi-coUinear singu- 
larities as 'q-singularities' for brevity. 

In Section ITT] we briefly recall DSM and describe our 
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extensions. Next, in Section III we address the problem 
of IR sensitive logarithms in the framework of the ACOT 
scheme. Finally, in Section [TV] we present a test of our 
method by an explicit numerical MC calculation of the 
charm structure function at NLO. 



II. DIPOLE SUBTRACTION METHOD WITH 
HEAVY QUARKS 

Let us consider a calculation of n-jet cross section to 
NLO accuracy (our notation is close to that of Ref. |11|). 
First, the LO contribution is (schematically) 

ai^O) ^UY^fa® j d^n,a |X„,a|' (1) 



2 



where fa is a PDF for parton a, the symbol (E) denotes 
convolution and Af is an appropriate normalization fac- 
tor. The n-particle phase space (PS) is denoted by d^n,a, 
while Mn,a is a tree-level matrix element (ME) with n 
partons in the final state. The jet function F„_a pro- 



duces the pertinent observable upon integration over the 
phase space. It has the property that i^„+i_a = Fn^a in 
the singular regions of phase space (IR safety condition) . 
Consider now the NLO contribution. Within DSM it 
reads 



(2) 



The notation A^ri°a°^'' corresponds to virtual correc- 
tions to |A^„,a| • The auxiliary function 2?„.a (called 
'dipole function' in the following) is chosen in such a 
way that: i) it exactly mimics soft and q-singularities 
of lA^n+i.a.r without doublc counting them in the over- 
lap region ii) it can be analytically integrated over 
the one-particle subspace dcfia defined symbolically by the 
relation d^n+i.a = d^n.a ^ d^a- The novelty here is that 
also the initial state q-singularities are considered when 
constructing the dipole functions. By construction, the 
first square bracket (real emissions) is integrable in four 
dimensions, thanks to the properties of jet functions. In 
order to leave the cross section unchanged we have to 
add the same dipole contribution as we have subtracted. 
This time, however, we use the integrated dipole function 
using its property ii). In general, this integrated dipole 
contains IR poles (in dimensional regularization) , which 
eventually get canceled against the poles in Mn°°^\ The 
resulting expression is, however, still not IR safe and the 
remaining q-singularities have to be factorized out. It is 
achieved by subtracting following 'counterterms' 



Cn.a=^Fab^\Mr 



(3) 



The quantities J-at are the renormalized densities of par- 
tons b inside parton a. We shall come back to these 
objects later in Section [Till 

Let us now define the dipole function more precisely. 
It is given by the sum of three (for DIS) classes of con- 
tributions 



n+1 

+ d\% (Pa; Pb) + ^S!k Pc) , 
fc=l ^ 



(4) 



where pa is the initial state parton momentum and pi 
with I = l,...n -|- 1 are final state momenta satisfy- 
ing momentum conservation q+ Pa = J2iPi- The three 
classes of dipoles are: (A) initial state emitter with fi- 
nal state spectator (lE-FS), (B) final state emitter with 
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FIG. 1: Kinematics for FE-IS (upper) and lE-FS (lower). 



initial state spectator (FE-IS), and (C) final state emit- 
ter with final state spectator (FE-FS). The 'emitter' and 
'spectator' partons are defined here following Ref. [81. 
The FE-FS case (C) does not involve initial states and 
was completely covered in [llj, thus we do not consider 
it further. 

The remaining dipoles in (E| read 



{Pa-Pi) -ml 



{M„ {pa^; Pa) \ V}^a^^,, \Mn {Pa^. P a) ), (5a) 



D\%{Pa]pB) = - 



{p,+p,f-ml X 



{Mn iPa; Pb) I ifl,,^ a\Mn {Pa, P b) ), (5b) 

where |A^„) is a vector in helicity and color space such 

that |A^„|^ — {A4n\-Mn) is a matrix element squared, 
summed/averaged over colors and spins (see [8J for de- 
tails). The kinematics is depicted in Fig. [l] For lE-FS 
the partons a, i are replaced by a single parton ai, while 
for FE-IS i, j are replaced by ij. Also the spectators ac- 
quire new momenta. All the new momenta in A^„ are 
constructed to be on shell and satisfy the momentum con- 
servation, also away from the IR limits. We adorn these 



3 



new momenta with tilde. The resulting sets of final state 
momenta in A^„ are denoted by pA and ps, respectively. 

The X variable is the Sudakov longitudinal fraction of 
Pa or pai along pa and the color operators T; pick up rele- 
vant color factors and account for color correlations. The 
dipole splitting functions V act in the helicity space and 
account for spin correlations. In the quasi-coUinear limit 
they tend to massive generalizations of splitting matrices, 
while in the soft limit they collapse to eikonal factors. 



A. Dipole splitting functions 

Here we present the realization of in a general mas- 
sive case. It is worth mentioning that our results extend 
those of Ref. [TT], and coincide in suitable limits with 
the partially massless cases considered there. 

First, let us define the tilded momenta as follows 



lE-FS: p'; = WAV''-iiAP^, p_ 
FE-IS: pt.^-^BV^-UBP'i, P^^n, 



Q^, (6a) 
Q^, (6b) 



where V = Pi +Pj , Q = V —pa (Fig- [l]) and u, w variables 
are calculated from the on-shell conditions. In the soft 
limit {t — > 0, — > 1. Moreover u — 1 — a; in the quasi- 
coUinear limit. 

Let us now list our dipole splitting functions starting 
with lE-FS case. For a gluon emission Q — > Qg (or 
Q ^ Qff) we have (with mq = ma — ruai, mi = 0) 



The correlation tensor V^V^ is defined by 



v^^ii-s)p';-sp^ 



+ [ml,-m'^+r'{l-2~z)]l^r^. (10) 

It fulfills trans versality relation pai ■ Va — 0. 

Now let us turn to the FE-IS case. Here the difference 
with [11 is that the initial state spectator can be massive. 
For Q — > Qg with = 0, mj — mij = mq we have 



2 

^^-32 ^ + 

ubVq + z 



(1 - e)z-2- 



Pi ■ Pj 



(11) 



Further for g — > QQ with m-ij = 0, = mj = mq we 
define 



V g-^QQ,aJ 



with 



V^ = zpf-(l-z)p^^ 



_ Usma , u. _ UN 



(12) 



(13) 



satisfying • Vb = 0. 

Finally for the g — !■ gg splitting we have 



UAV] + Z 



-^+{\-e)uA-2- 



1-MA 'Tig 
~Z Va 



(7) 



where z = pi ■ Pa/Va with Va = V ■ Pa- Further, is a 
mass scale needed in Z? = 4 — 2e dimensions, and Vj = 



1 — {mamj /^y)^ with j ~ pj ■ pa- Above (and in what 
follows) the unit matrix in helicity space is suppressed. 
For the g QQ splitting we get 



y(A) 



1 - 



1-e 



2uA (1 - ua) - 



l-UA 

~Z Va 



(8) 



with the configuration of masses ma = 0, mq = m^i — 
mi. Note that although the initial state parton is mass- 
less, this splitting is not present in [111 IT2j . 

Finally for Q — )• gQ, where a heavy quark is radiated, 
we have (with m^ — mi — mq, mai = 0) 



'Q^gQ.jj = 87r^?^a,CF (1 - e) 
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-g''^ (1 - ua) + 



1 - UA 22Va - 



Q 



(9) 



I - z + ub 



z + ub 



(14) 



Despite the massless splitting the above function gains 
some massive factors from spectator via the vectors. 



B. Dipole integration 

The phase space factorization required in Eq. ^ is 
actually realized as a convolution in u treated as a free 
variable. In this context we denote it by u (without 
tilde). The construction requires dropping on-shell con- 
ditions for Pj in lE-FS or pij in FE-IS cases and fixing 
two invariants X , y, such that we can determine w = 
w (m, A", 3^). Let us denote the corresponding off-shell 
vectors as p^ (u, X, y) = pj (u), (m, X, y) = p^^ {u). 

Our version of PS factorization formula has the form 
(for lE-FS case) 



jdud<^,, {Q,pM:Pa) dc^i%^^^ (u) , (15) 



4 



where the arguments of the reduced phase space 

are the incoming and outgoing momenta (separated by 

semicolon). The subspace measure is 



, (16) 



4(27r) 



where v = -y^l — ra^V'^ /V^ and dQD-2 is a solid angle 
on transverse hyperplane in the Q-Pa CM system. The 
jacobian reads 



J 



du 



(17) 



with 
fixed 



subscripts corresponding to variables kept 
during differentiation. The bounds on 

± 2vp) , where 

and the barred 



1 - 2 I 



i 1 



(?7if - ml 



masses are rescaled hy . 

An analogous formula can be obtained for the FE-IS 
case. One has to replace the set by ps and make 
replacements: pai ^ Pa in d^^n and pj — !■ pij in the 
jacobian . 

The choice of invariants X , y affects the PS generation 
in (i$„ and 'plus-distributions' handling in the integrals 
over d4). In the following we take X — -^^y — Va- 

We note that in the desired integral J dcfiD, with D 
given in only V and the propagators can depend on 
z. Further, it can be shown that, thanks to the gauge in- 
variance and properties of ^ , the helicity correlations 
in (|5| vanish after integration. Therefore the integral 
over dr2£)_2 is trivial and we are left with a color cor- 
related matrix element and a scalar function, which e.g. 
for lE-FS reads 



(u) 



(A) 



a—>-att. J 



(Pa-Pi) --ml 



(18) 



where (. . .) denotes average over helicities in D dimen- 
sions. The formula for FE-IS is analogous. 

All the integrals have been calculated analytically in 
D dimensions [7J , and the details will be given in a forth- 
coming paper. 

As an example, let us present the result for the lE-FS 
3 — QQ dipole integral as we shall make use of it in 
Section |III[ In the D — A limit we get 



/(^) - (u) 



9 



Vq - + Vv2 (1 + 2p) 
8Tr»?q'?p2 (1 - u)p 1 



(19) 



where Pgq (u) 



Tr [1 — 2w (1 — u)] is the standard split- 
ting function, — to|/27, 7/q = mq/27 and other 
ri\ = X/27 for X = Va,V'^,J- Note that in the above 
formula a non-trivial dependence on u is hidden in p and 



III. FACTORIZATION OF 
(QUASI-) SINGULARITIES 

In general, the integrated dipole functions contain IR 
poles, which get canceled upon including virtual correc- 
tions. There are no further singularities in the massive 
case. Nevertheless, some IR sensitive terms (q-singular- 
ities) still remain in the lE-FS dipoles, and we factorize 
them out by subtracting counterterms (|3| according to 
the ACOT scheme, as we explain below. 

Let us consider the integrated lE-FS dipole in the limit 
of vanishing mass of the emitter or emitted parton (the 
spectator j can remain massive). For instance, using Eq. 
( 19 1, we obtain 



g^QQ,j^> 27r 



Pgq {u) log 



log 77q 



2rR,u(l -u) 



oiv^a), (20) 



The counterterms Cn,a are given in terms of the par- 
tonic densities J-abi which are defined via matrix elements 
of bilocal operators on the light-cone and can be calcu- 
lated perturbatively. The emerging UV singularities have 
to be renormalized which, in turn, specifies the evolution 
kernels for PDFs (see [TJj for a review) . 

According to the ACOT scheme, this standard proce- 
dure can also be used with massive quarks [6J. In par- 
ticular, one can still define partonic PDFs and renormal- 
ize them using MS, assuring standard DGLAP evolution 
equations (actually one uses the composite CWZ renor- 
malization scheme [TS] ) . For asymptotically large , the 
above construction leads to the standard MS calculation. 
On the other hand, when ~ mq the Cn.a countert- 
erms approximately cancel the a^^'^^ contribution, basi- 
cally resulting in the 'fixed order' description ^2^. 

Using the Feynman rules for PDFs [HI [TB] one can 
show that in the MS scheme 



27r 



Tr (l-2z(l-^)) log 



(21) 



2n 



Cp 



l + (l-^f 



log — 



2 log z - 1 



(22) 
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QQ 



^gg (^) 



27r 

1 + 
1 - z 



2^ 



cf 



21og(l - z)-l 



(23) 



2CA 



+ <5(l-z) 



1 



l + z(l-z) 



(24) 



where z is the fractional momentum of parton 6. 

We have checked that using Cn,a with the above Tab 
(intermediate steps require usage of color conservation 
for Tj operators) we get IR safe dipoles. Moreover, in 
the considered limit we get the same dipoles as those of 
Ref. |TT] with the standard massless collinear subtraction 
term (i.e. with Tab {z) = -UsPab {z) /27re). 

There is one comment in order. Certain ambiguities [6_ 
of the ACOT scheme stimulated construction of the so 
called Simplified ACOT (S-ACOT) scheme g], where all 
initial state masses are set to zero. This scheme is conve- 
nient for higher order calculations |17j, which otherwise 
are very cumbersome. Our approach is however differ- 
ent as it operates on more exclusive quantities, hence 
we resolve the ambiguities differently. Namely, we set 
the initial state masses to zero in cr^'"*-'' and in the Cn,a 
contribution to cr^^'^'-'). This results in massless integra- 
tion limits in the corresponding convolutions, in analogy 
to the standard DGLAP evolution equations. This ap- 
proach is not strictly required by our formalism, but pro- 
vides a better transition to the 'fixed order' description, 
cf. solid and dotted lines in Fig. |2] We point out that 
we do not set any other masses to zero, neither for initial 
nor for the final states. If, however, one wants to use the 
S-ACOT scheme for jets, then only the dipoles for the 
g — >■ QQ initial state splitting are needed. 
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FIG. 2: The charm structure function, F2c, calculated in dif- 
ferent schemes: GM from our MC code at NLO (red solid 
hne), Zero-Mass at NLO (black dashed line), Fixed Order 
(BGF) at 0(05) (blue dotted line). The calculations were 
done using CTEQ5 PDFs at xb = 0.05. 
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— , d[^X r\ ~ d'^^J- — as well as their integrals. Since 
g,Q,Q 9 s,Q,Q 

there are two lE-FS dipoles we need two collinear sub- 
traction terms with corresponding TgQ and ^qq- 

We observe, that the soft poles coming from virtual 
corrections cancel against the sum of soft poles coming 
from lE-FS and FE-IS dipoles. The result of numerical 
integration perfectly agrees with the semi-analytical for- 
mula calculated independently according to [21 [19]. In 
Fig. [2] we show that our result indeed interpolates be- 
tween the massless NLO calculation at high and the 
fixed order, 0{as), calculation at low Q^. Let us stress 
that this result is obtained by the MC integration of a 
fully differential cross section (2). 

This is clearly the first simplest test. It verifies only 
a part of the dipoles and collinear counterterms. More 
tests are due in course of the code development. 



SUMMARY 



IV. NUMERICAL TESTS 

We have partially implemented our method in a ded- 
icated C++ program. The MC integration and event 
generation is accomplished using the FOAM engine |18| . 

In order to test our approach, we have used our 
program to numerically calculate a quantity which can 
be obtained independently, namely the charm structure 
function i^2c- 

The LO contribution comes from 7Q — > Q (plus the 
same with Q). The real NLO corrections are 7Q — > Qg 
and 75 — >■ QQ, while the virtual ones are given e.g. in 

m- 

In order to perform the numerical integration we need 
three distinct dipoles D^^^ = D^^^ 0^%^^ = 



In the present paper we have briefiy described — with- 
out going into technical details — a general mass scheme 
for jet production in DIS. The scheme is based on the 
ACOT factorization theorem which is proved to all or- 
ders. The proof is constructed for structure functions, 
but it actually holds for all IR safe quantities. This enti- 
tles us to use this factorization scheme for jets, provided 
we have a suitable dipole subtraction method that can 
deal with heavy quarks in the initial state. We have con- 
structed such a method, generalizing the existing results. 
We have checked that this method removes all potential 
collinear singularities as desired, and that the remaining 
dipoles coincide with those of the massless calculation in 
MS. 

All the crucial calculations are already performed j7J 
and the details shall be presented in a separate publica- 
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tion. The computer code suitable for the jet production 
in DIS is under construction. 

An extension to hadron-hadron collisions is, in princi- 
ple, straightforward. Formally, one only needs to intro- 
duce an additional class of dipoles for the initial state 
emitter and initial state spectator. Nonetheless, we are 
aware of challenges arising at higher orders (above NLO) 
— like the feasibility of dipole subtraction method or 
non-cancellation of soft singularities when two initial 
state partons are massive. These issues have been widely 



discussed in literature, but practical procedures still have 
to be analyzed. 
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